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G\ 

Q\ • Electron transmission and wave functions through junctions with a pair of a pentagonal defect 

0^ | and a heptagonal defect connecting two metallic carbon nanotubes are analyzed by the analytical 

calculation with the effective mass equation. The energy region \E\ < E c is considered where the 
5^ ' channel number is kept to two. Close relation between the transmission rate and the wave function is 

found; the transmission rate is given by the inverse squared absolute value of the wave function. The 
dependence of the transmission rates on the energy and on the size of the junction is clearly explained 
by the nature of the wave function. Though the wave function and the transmission rate calculated 
by the tight binding model agree well with the corresponding analytical results by the effective 
mass approximation, the discrepancy becomes considerable when \E\ ~ E c . To study the origin of 
this discrepancy, an efficient numerical calculation method is developed with a generalized transfer 
matrix for the tight binding model. Their numerical results are compared with the corresponding 
analytical ones and the results show that the origin of the discrepancy comes from the evanescent 
, waves with the longest decay length in the tube parts. 
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I. INTRODUCTION 



Recent experimental development on the carbon nanotubes, especially electronic transport measurements for indk 
Ch ' vidual nanotubesjj has promoted much interest in the nanotubes as one-dimensional conductor with nanometer size.EI 
Many experimental works concern the theoretical prediction that the nanotube becomes metallic or semiconducting 
i ^ i ■ accordiiag-jto its circumference.Eru The conductance of the metallic nanotubes with potential energy disorder has been 
studied.M The junction structures with a shape p£ a part of a cone connecting two nanotubes with different radii have 
. been akp|Observedj3 and studied theoreticallyETtll They are formed by a pair of a pentagonal defect and a heptagonal 
defect Jl3~llil By composing the nanotube junctions, the electronic circuits with nanometer size might be designed. The 
atom bond network of the nanotube junctions is uniquely determined by its development map with the vector of the 
' circumference of the connected thicker nanotube i?5 and that of the connected thinner tube R7. We have calculated 
(f-) , the conductance with Landauer's formula, i.e., the transmission rate, for the various junctions connecting two metallic 
nanotubes by the tight binding modeLBs Surprisingly, the dependence of the transmission rate on the parameters, 
G\\ R 5 , R 7 and E, has a very simple form obeying the scaling law; in the energy region \E\ < E c where the channel 
number is kept to two in both the tubes, the transmission rate is independent of detailed atomic arrangements as well 
as the angle between the two tube axes but determined only by the two parameters; one of the two parameters is the 
ratio of the circumference of the tubes R7/R5 and the other is the scaled energy E/E c . The close relation_b.etween 
the transmission rate and the wave function in the junction part is also found by the tight binding model.mil When 
. E = 0, the transmission rate decays with the thickness of the junction as 4/(2 + (R 5 /Ry) 3 + (R7 / R 5 ) 3 ) while the 
spatial decay of the corresponding wave function in the junction part obeys the same power lawH When R7/R5 ~ 0.5, 
the transmission rate shows a peak structure as a function of the energy, while the corresponding wave function shows 
a resonant feature, i.e., its amplitude is enhanced in the middle of the junction partjj _ . 

Recently, Matsumura and Ando have confirmed the power law decay for E = by using the effective mass theory.E£l 
We generalized their discussion to more general energy region \E\ <rEc, and obtained the complete analytical form 
of the transmission rate with the two parameters, R7/R5 and |^|/J5 c Ji3 The agreement of the analytical transmission 
rate with the numerical one is found to be fairly good so far as \E\ is not very close to E c . By our generalization of 
the effective mass approximation, the .band structures of the periodic multiple nanotube junctions, which are called 
the 'helically coiled nanotube' in Ref.tll , can be also obtained as closed analytical forms and their .relation to the 
configuration of the pentagon and the heptagon are clearly explained based on the symmetry features. E£l The effective 
mass theory is derived from the tight binding model and is valid when the energy is near zero. It is suitable to 
obtain the analytical results and to explain the origin of the scaling law because it treats not the discrete..sp.ace but 
the continuous space, i.e., it does not necessitate the detailed atomic structures of the honeycomb lattice.BEjll3'Ea~E2l 
Though the analytical transmission rate of the nanotube junction by the effective mass theory has been studied in 
detail in this way, the close relation between the transmission rate and the wave function observed in the numerical 



1 



result is not clarified by the effective mass theory yet. In this paper, the spatial variations of the wave function is 
expressed unambiguously by the effective mass theory and it gives an intuitive explanation for the dependence of the 
transmission rate on Rt/R 5 . Furthermore it will be also found that its dependence on the energy has close relation 
to the other parameter \E\/E C of the scaling law. These points are discussed in section III. The concise expression 
of the wave function characterized by only the two parameters will be useful when the STM images of the nanotube 
junction are analyzed. 

On the other hand, the discrepancy between the numerical transmission rate with the tight binding theory and the 
analytical one by the effective mass approximation become considerable when |J5| is very close to E c ; the numerical 
transmission rate shows a sharp dip there while the analytical one does not show. Since the effective mass theory is 
an approximation derived from the tight binding model, the numerical result by the tight binding model is considered 
to be correct one. To discuss the origin of the sharp dip, the wave functions calculated by the tight binding model 
have to be compared with those calculated by the effective mass theory. To calculate the wave function by the tight 
binding model more efficiently, a new numerical method is developed in section II. From the comparison of the wave 
function calculated by the method in section II with the analytical one, we found in section IV that the evanescent 
waves in the tube parts and the defect levels at the pentagon and the heptagon have important roles to form the dip 
structures. 



II. THE CONDITIONED TRANSFER MATRIX METHOD 

Fig|l] shows the development map of the nanotube junction. It is characterized by bars representing zig-zag segments 
of the C-C bond network in the circumferential direction of the single junction. They are aligned and numbered along 
the direction ei — e*2 where ei and e*2 are the basic translation vectors of the graphite plane. Each bar is connected 
with adjacent bars by the remaining C-C bonds. The network is rolled up so that the atoms denoted by j in the 
bottom of the j'th bar and the top one of the same bar denoted by j' become the identical one. The j and j' sites are 
shown by the filled and the open circles in Fig.pl respectively. When j < 0, the j'th bar is formed by the (to + n) unit 
cells of the 2D graphite, where the m unit cells are aligned along ei and the others are aligned along e 2 . So the (to, n) 
tube defined in Ref.Q is formed there. When 1 < j < h, the j'th bar is defined as such that made by adding (j — 1) 
unit cells and one atoms along the direction of el to the bottom of the zeroth bar. For j > h + 1, j'th bar is made 
from the Zi'th bar by exchanging (Z2 — 1) unit cells and one atoms in its bottom aligned along el with I2 unit cells 
aligned along e%. The network made in this way represents the junction of the (m,n) tube (j < 0) and the (rn2,ri2) 
tube (j > li + 1) , where l\ = bpn% + n 2 ) — (to + n) > and l 2 = n 2 — n > 0. This junction is called an (to, n)-(m2, n%) 
junction hereafter. Then Fig.|l| corresponds to the (2,2)-(2,5) junction. A 7-membered ring is formed at the bottom 
of the 1st bar and a 5-membered ring is introduced between the Zi'th bar and the (l\ + l)'th bar. There are only 6 
mcmbercd rings elsewhere. 

The number of the bonds connecting the (j — l)'th bar and the j'th bar is denoted by bj as 

bj = to + n (j < 0) 

bj=m + n + j — 1 (1 < j ' < li) 

bj = m 2 + n 2 (h + 1 < j) • (1) 

The amplitudes of the wave function in the j'th bar are represented by the vector Cj = (Cj,i, Cj t %, ■ ■ •). Its component 
Cj t i can be classified to the two groups according to whether the corresponding site connect with the right (j + l)'th 
bar or with the left (j — I)'th bar. These two groups are denoted by Tj and lj, respectively, as is shown in FigJ|. Since 
each site in j'th bar belongs to only and necessarily one of the two groups, the dimension of Cj equals bj + bj+i, which 
are denoted by dj hereafter. The bonds between the (j — I)'th bar and the j'th bar are numbered as FigJ| from the 
bottom to the top. According to this way of numbering, the amplitudes are also numbered by rj-i.i and lj^, which 
have the bond numbered by i (i = 1, 2, — , bj). Then the tight binding equation become 

Erj- hi = j(lj,i + (2) 

i' 

and 

El jti = j(rj-i,i + Y c hi") ( 3 ) 

i" 

where the first and the other terms in the right hand sides represent the bonds between the neighboring bars and those 
within the bar, respectively. Here we use the tight binding model including only 7r orbitals with common hopping 
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integral, 7(~ — 2.7eV), and common site energy chosen to be zero. In order to get more accurate results, the effects 
from the mixing with a orbitals caused by curvature of the graphitic plane also have to be considered. But we believe 
that this tight binding model gives a transparent view about this system and qualitative valid results. Another reason 
why this model is used is that the purpose of this paper is to focus on the effects from the connectivity of the bond 
network rather than those from the curvature. 

The tight binding equation, eqs. @) and (^) can be summarized by the matrix form as 



B 3 c 3 =Q 



(4) 



where Aj and Bj are a 2b j x dj matrix and a 2b j x dj+i matrix, respectively. In the left (right) tube pars, i.e., when 
j < (j >l% + 2), these matrixes become constant matrixes, Al and Bl (Ah and Br). The transfer matrixes for the 
tube parts are obtained as Tl = —A~£ Bl for the left tube and Tr = — A~^~Br for the right tube. The eigen values 
and eigen vectors of T M are classified into two groups as {/3f,-0+i} an d {l//3f,^j}. Hereafter, \i — R and fj, = L 
represent the left thinner tube and the right thicker tube, respectively. The propagating waves in the former group, 
which are tp+i with |/3j| = 1, are assigned to i = 1, • • • , n^, where n M is called the channel number of the corresponding 
tube. They carry the probability flow with the positive velocity, m. The propagating waves in the latter group, 
ip-i — ip* , carry that with the negative velocity, — «j. From now on, the normalized extended states with unit flow, 
i>i = ipil 1 \fv~i an d ip—i = ip-i/y/vi are used instead of tpi and ip. 

and the sign of i means the direction along which they decay exponentially, i.e., \/3jf\ < 1 

The wave function in the tube parts are represented by superposition of these eigen vectors and its coefficients are 
represented by x as 



The other states for \i\ > are evanescent waves 



(5) 



where the z'th column of U± is the normalized eigen vector and A M is the diagonal matrix whose (i, i) element is 
the eigen value /3f . 

Unlike the tube parts, the 'usual' square transfer matrix cannot be obtained in the junction part, because the 



dimension of Cj and that of c +\ are different, i.e., the Aj 
inverse matrix. Nevertheless there is a dj 
by Aj and defined by 



for 1 < j < l\ becomes a rectangle matrix which has no 

It is denoted 



x 2bj matrix which is similar to the inverse matrix of Aj 



Aj = (A\AjY 



(6) 



This matrix satisfies the condition AjAj = 1 but AjAj ^ 1. By this pseudo-inverse matrix, we can define the 'transfer 



matrix', To 



-AjBj. If eq.(W) is satisfied, then 



Cj-i — TjCj 



(7) 



is satisfied, but its converse does not hold generally. It means that Cj+i in eq.(Q) cannot be chosen to be arbitrary 
but has to satisfy 



(S J +^T J )c J =0 



(8) 



To make eq.(^|) equivalent to eq.(^), the condition (^) is necessary, so we call Tj a 'conditioned transfer matrix' 
hereafter. In our problem, the number of the independent rows of (Bj + AjTj) is only one, so cq.(||) can be written as 



: s 3 ■ Cj = 



(9) 



where * Sj is one of the nonzero rows of (Bj + AjTj). Multiplying Tj generates the coefficient at the bar with the 
decreased j, so we define the positive direction along which j decreases, i.e., from the right thicker tube to the left 
thinner tube. 

By using the inverse matrix 4 (' V^) = (Z7+, U 1 !)^ 1 and Kj = TjTj + \ ■ ■ ■ T/ 1 , the relation between and 
the right going waves in the thinner tube ar^ are summarized as 



'4+1 



V 4 s 2 K 2 J 

















C/i + l = 






(10) 



V o J 
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The (rri2 + n%) x (2m2 + 2ri2) matrix in the left hand side of eq.([lO|) is represented by Q hereafter. 

Now we define all the necessary things to calculate the transmission rate and the wave function. When the electron 
is incident from the right and transmitted into the thinner tube , i.e., when x R ^ and x^L = 0, the reflected wave 



x R and the transmitted wave arr are obtained from 



x R = -{QU R )-\QU R )x R (11) 

and 

x L + = V^K a (U R - U R {QU R )-\QU R )) x R , (12) 

respectively. The corresponding wave function in the junction part, i.e., Cj for 1 < j ' < l\ is given by 

Sj = Kj (U R - U R {QU R )- 1 {QU R )) x R . (13) 

For the inverse direction of the incident electron, x R = and xb. ^ 0, the reflected wave x+ and the transmitted 
wave x_ are obtained from 

x% = VfKoU^QU*)- 1 ( ^ ) , (14) 

and 



x> R = (QU R )- 1 ^} . (15) 

In this case, the wave function in the junction part is given by 

S i = K j U?(QU?)- l ^ y ) . (16) 

The transmission rate is obtained from eq. (|l^) or eq. (|l^) , while the reflection rate is obtained from eq. ( |j"l| ) or eq. (|l4| ) 
from the matrix elements corresponding to the open channel, i.e., x^ with i = 1 ~ n^. We have confirmed that 
unitarity holds very well in the conditioned transfer matrix method. Agreement of their results with those calculated 
by usual recursive Green's function methods is also quite good. In both the methods, matrix inversions are necessary 
about \\ times, but they can be done much faster in the conditioned transfer matrix method than in the recursive 
Green's function method, because matrixes which have to be inverted are real sparse symmetric matrixes, A^Aj, in the 
former method while those in the latter method are complex matrixes with no symmetry. Therefore the 'conditioned 
transfer matrix' method is much faster than the recursive Green's function method. This efficiency is expected to be 
more important when more general related problems are calculated, e.g., those including electron-phonon interaction 
or electron-electron interaction self-consistently. Another advantage of the conditioned transfer matrix method over 
the recursive Green's function method, which is not essential but practically important, is that it is more intuitive 
and easy to be implemented. 

In section IV, the wave functions are calculated numerically by the conditioned transfer matrix method and com- 
pared to those calculated analytically by the effective mass equation. 



III. EFFECTIVE MASS THEORY AND ITS APPLICATION TO THE SINGLE NANOTUBE JUNCTION 

Fig. U shows a development map of the nanotube junction neglecting the atomic sites. The vectors R§ and 
i?7 represent the circumferences of the thicker tube and the thinner tube, respectively. Lines -E7P7P5-E5 and lines 
F7Q7Q5F5 are sticked each other so that the points connected by Rj become identical (j — 5,7). Then Ps(= Q5) 
and P-j{— Q7) turn out to be the centers of a heptagonal defect and a pentagonal defect, respectively. These points 
are marked also in Fig.[l]. Thus from now on, the indices '7' and '5' are used to represent the thinner and the thicker 
tube, respectively. 

The equilateral triangles 'AOP7Q7' and 'AOP5Q5' with bases-y-VC^' and 'P5Q5' have a common apex O, which 
is chosen to be the origin of the coordinate (x,y) in this paper .113 The origin O lies between — (m + n)'th bar and 
— (m + n — l)'th bar , since Q7O = me*2 + n(ei — e*i). 
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Then the position of a general site, f, can be labeled by (3,(71,(72), where qi and q 2 are integer components 
representing the position of the unit cell, q = qi ei + q 2 &2, and s represents the sublattice, s = A, B. Relation between 
the position vector of the site, f, and its label (s, qi,q 2 ) is written by 



r%A, qi,qz) = q+ (ei + e 2 )/3 
\/3, 

T (?l + 92 + - 3/ -. . 2 



V^3 2 1 

X (?l + <22 + -z)&x + ^{q2 - qi)e y (17) 



and 



r(B,qi,q 2 )=q + 2(ei + e 2 )/3 



= +92 + ^)e x + i(<7 2 - ?i)e y . (18) 

In the above e^ = (ei + e*2)/\/3 and e^ = e 2 — ei, which define the coordinate system (x, y) shown in Fig.[l]. The 
amplitude of the wave function at these sites is denoted by ip(s, qi, q 2 ). In the effective mass theory, the wave function 
is represented by 

ip(s,qi,q 2 ) =F? {f)w iq '- q ^ + Ff {f)w^- q ^ {s=A,B). (19) 

Here w = exp(z27r/3) and r is defined by eqs.(|l7j) and (18) as a function of the label (5,(71,(72). In eq.(19), F^'g , 
^(91-92) anc j ^(92-91) are envelop wave functions, the wave function of the Bloch states at the K and the K' point 
, respectively. This wave function is expressed by a vector, ip = {F£jjL\—Fg(r J r e x /V3),F^ (r),Fg (r + e x /v3)), 
hereafter. This definition of F's is different from other referencesBOEjlij by certain factors. The reason why this 
definition is used is that the representation of the time reversal operation /, lip = ip* , becomes simpler as 

I(FZ,F*,F? ,*#') = «Ff )*,(F*' r,(F«)*,(FK)*) . (20) 

By using eq. ( |l9|) , the effective mass equation are derived from the tight binding equation used in the preceding section. 
When energy E is zero, w qi ~ q2 and w q2 ~ qi are solutions of the tight binding model, so that the corresponding solutions 
of the envelop functions F are constant. Thus when E is not zero but close to zero, spatial variation of the envelop 
functions is slow compared to the lattice constant, | et,. | = \e v \ = a ^ 0.25nm. In this case, it is a good approximation 
to take only the first order term-m-tiLC Taylor expansion of the envelop function as F(r + e x ) ~ (1 + ad x )F(r). From 
this approximation, one obtainsQotJ 

(-id v + d x )F*(r + b) = eF« (r) (21) 
(-id y - d x )F«(r) = eF«(f+ b) . (22) 
{id y + d x )F«'{r + b) = ePf (f) (23) 

(id y -d x )F¥(r) = eF*'(r + b) , (24) 

where e = 2E / (\/3^a) and b = e x /y/3. When the plane wave solution Fj£ , Fg oc exp(±zfc • r) is used in eqs.(^i"|) and 
(p2|), one can get linear isotropic dispersion relation, 

k=\k\ = \e\ , (25) 

and Fg(r + b)/F^(r) = ±exp(z?7), where (k x ,k v ) — (— e sin 77, e cos 77), i.e., 77 is the angle of k with respect to the y 
axis measured anti-clockwise. It follows that the corresponding wave function, ipK±i is written as 

$K± = (e-"' /2 ,±e"' /2 ,0,0)e ±l ^ . (26) 
The wave number must satisfy the boundary condition in the metallic nanotube parts as 

k-R 3 =2^l, (j=5,7), (27) 
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where lj is an integer representing the number of oscillations of the envelop function around the circumferences. As for 
the semiconducting tubes, the boundary condition becomes different from eq.(p7|)Ej, but we concentrate our discussion 
to the metallic nanotubes in this paper. The upper sign and the lower sign in eq.(p6|) represent the direction of the 
propagating waves. When E is close to zero, i.e., k ~ 0, only possible number of lj is zero. It means that k is 
perpendicular to Rj so that 77 is the angle of Rj with respect to x axis. Within the effective mass theory scheme, the 
possible maximum value of I5 (£7) is the channel number of the thicker tube (the thinner tube). Therefore the range 
of the energy where the channel number is kept to two in the thicker tube ( the thinner tube) is I-E/7I < \J~Zna/ R§ 
(I-E/7I < y/3na/Rf). Discussions in this paper are concentrated to the energy region where the channel number is 
kept to two in both the tubes. It is represented by \E\ < (\/37r|7|)(a/i?5) = E c and E c is called a threshold energy 
hereafter. From the propagating waves near the K point, ipK±, the other propagating waves ipK'± are obtained by 
the time reversal operation ( po|) as 

$ K ,± = lj KT = (0, 0, e">/ 2 , Te -«,/2 )e ±*(k-f) (28) 

Note that the direction of the propagation is reversed by the time reversal operation /. In order to discuss the wave 
function in the junction part, the polar coordinate (r, 6) is useful. Its relation to the coordinate (x,y) is the usual 
one, i.e., r = \J x 2 + y 2 , tan 6* = y/x. Then the wave function satisfies the wave equation r 2 (d 2 + d 2 + k 2 )F — 
{z 2 d 2 + zd z + dg + z 2 )F = 0, where z — kr. The solution is represented by Bessel functions J m and Neumann 
functions N m as 

00 

F= J- e m9 (c m J lml (z) + d m N lml (z)) . (29) 



m— — 00 



In Fig.Q, site i is identical with site i' for i = 1, • • • / 1; while corresponding labels (s, qi, 52) are different between i 
and i'; that for i is (£>, 1, —4 — i) and that for i' is (A, 3 + i, —2 — i). In this way, the wave function in the junction 
part must satisfy the condition 

i;(A,q 1 ,q 2 )=rp(B,q 1 + q 2 ,-q 1 -l) . (30) 

Here transformation of the label (A, qi, q 2 ) — * (B, qi+q2 + 1, —91 — 1) is equivalent to clock-wise rotation by 7r/3 with 
respect to the origin O. From eqs. ( |30| ) and (|l9|), the boundary conditions in the junction part are derived as 



and 



The other boundary conditions 



and 



(r,9 + Tr/3) = wFK(r,9) , (31) 



Ff(r,0 + 7r/3) = -F#(r,0) . (32) 
w 



Fg(r,e + Tr/S)=wF^(r,e) , (33) 



F*'(r,0+7r/3) = -Ff(r,0) . (34) 
w 

are also obtained in the same way. The same boundary conditions are firstly discussed by Matsumura and Ando.EH 
Difference of eqs. (|3l]) - ( |34| ) from those of Matsumura and Ando by certain factors is due to the difference of 

the definition of F A ' B . From eqs. ( |3~l| ) and (^), terms in eq. (^9|) for F% and Fg are not zero only when 
m = 3p + 2 (p =integer ). Because the open channel L — in the tube parts is spatially uniform along the 
circumference, it is better fitted to the components with smaller \m\ in eq. ( |29| ) than to those with larger \m\. So we 
assume that one can neglect all the terms except those with p = and p = — 1 in eq. (|2^) (Assumption I). Then the 
wave functions can be written as 



Ff =e M / 2 (z) + e-' fl /i(4 (35) 



and 



G 



n 2i6 



Hz) 



'AM 



(36) 



where 



f m (z) = c m J m (z) + d m N m (z) (m = l,2). 
From eqs. ( ^l| ) and fl24|), the other two wave functions i 7 ^ and F^" can be derived from F4' and Fg as 

(-e l9 / 2 (z) + e- j2e /iW 



^ s "Id 



(e ie f 2 (z) + e- i2e h(z)) , 



(37) 

(38) 
(39) 



where 



fi(z) = Cl J 2 {z) + d x N 2 {z), 

Hz) = c 2 Ji(z) + d 2 N x (z) . (40) 

In the above we used the recursion formula of the Bessel functions and Neumann functions.0 It is easily confirmed 
that eqs. ( |38| ) and (|39|) satisfy the boundary conditions eqs. (^) and (|34|). The amplitude of the open channel in the 
tube, which is denoted by a, is obtained from eq.(|26|) as 



2,/Ri 



dx {j) {e 1 ^ F% ± e^ff) (j = 5,7) 



(41) 



for the K point. The value of \a\ 2 equals the probability flow as shown in Appendix. The indices + and — mean the 
direction along which the electron waves propagate. The path of the integral of eq.(^l|) is the straight line PjQj, the 
angle of which with respect to x axis is denoted by r/j. Equations for aj± are obtained from eq. (^) by replacing ±, 
rjj and K in the r. h. s with =p, —r/j and K' , respectively. To simplify the calculation, the integrations in the above 
are transformed as 



dx u) -» Ra 



-f+Vj 



(42) 



If variation of the wave function along the radial directions is slow near r = Rj, this replacement can be allowed 
(Assumption II). The relation between the amplitudes of the open channel in each tube, dy = * (aj+, c^j+j a f-> a f-)i 
and the coefficients representing the wave functions in the junction part, g = 1 (c 2 , d 2 , c\, d\), are summarized in the 
followings. 



= -VRlPivjWA^LikRi)^ Y 3 g , 



(43) 



where M is a constant matrix given b 

M = 

A(rf) is a diagonal matrix, where Au = A3 

P(V) = 



( 







V3 
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(44) 



P{rj) is defined by eq. as 



(45) 
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The matrix elements of L(z) are Ln = L33 = Ji{z), = L34 = N\{z), L 2 \ = £43 = ^2(2) and L22 = £44 = N 2 {z). 
The other matrix elements of L(z) are zero. From eq. (|43|), the relation between CS7 and d?5 is given by a 7 — Y 7 Y^~ as 
with the three parameters, kR 7 , kR§ and </> ee 777 — 775 (angle between i?5 and i?7 in the development map). This reads 
as 



"7+ \ _ ( h, t 2 \ ( a 5+ 
a 7 - J \t 2 , tl J \ a 5 _ 



where 



(46) 



tl ~ h +\ iBin(l0), cos(fV) I ' (47) 



* 2 = *- ( ;rftt iSSf' i ■ (48) 



isin(|</>), cos( 2 
The factors h+ and in eqs.(^) and ( f48| ) are represented by 

h± =-i(Xi 2 (fci?7,fci? 5 )T^2i(fci?7,A:i?5)) + ^= (xix{kR 7 ,kRs) ±^X 22 {kR T ,kR 5 )^ (49) 

where 

Xij^i, z 2 ) ee ^/m^J^i)^^) - N l {z 1 )J J {z 2 )} . (50) 

The considered energy region \E\ < E c corresponds to the region, < kR$ < 2tc. It can be easily confirmed 
analytically that eq. ( [46"| ) satisfies the time reversal symmetry and unitarity. The transmission rate per channel 
denoted by T is calculated from eq.([l6|) as 

T=l/\h+\ 2 = = = , (51) 

{(1/6)E? =1 Eh(3/4F- 2 4(^,^)} + 1 

and it gives the conductance a as a — 2T by Landauer's formula. In Fig. [|, the solid lines represent the values of T 
calculated by the tight binding model, and the dotted, line and the dashed line represent those calculated by eq. j5l]). 
Agreement between the two methods is fairly good£3 When E = 0, T becomes 

4 

T = (R 5 /R 7 r + (R 7 /R 5 r + 2 • (52) 

Eq. (|5^ ) reproduces well the numerical results in Reference. 1. 

We are now ready to consider the analytical calculation of the wave functions by the effective mass theory and its 
relation to the transmission rate. Let us consider the case where the two waves cfg? (i = 1, 2) are incident from the 
thicker tube, and there is no incidence from the thinner tube, i.e., a 7 ~ — 0. As in section II, it is called the positive 
incidence direction. Whenever the two incident waves are 'orthogonal' and have the same flow, ii.e., \a£l\ = \a^_\ 
and d?5_|i • d? 5 + = 0, the corresponding two transmitted waves are also "orthogonal" and have the same flow, because 
t\ti — 1/T. Note that it does not hold generally and characteristic of the nanotube junction. 

The two orthogonal transmitted waves are represented by ay~2 =* (ai,a 2 ) and effl, =' (— a^al). In the analytical 
calculation, we use different normalization from that of section II; we take the transmitted waves with unit flow 

represented by |S 7 l ?| = 1 in this section while the incident waves have unit flows, |<3 5 +| = 1, in section II. Corresponding 
wave functions denoted by ip^ are obtained from g = Y^~ l a^ and eqs.(35)-([fb"|). They depend on the phases of a\ 
and a 2 , but the sum of the squared absolute values of them, ^(s, qi, q 2 ) ee |^>W(s, qi, q 2 )\ 2 + |t/^ 2 '(s, qi, q 2 )\ 2 , does 
not depend on them and can be determined uniquely. Therefore we discuss "J = l^ 1 ^ 2 + |?/>( 2 )| 2 rather than l^ 1 )) 2 
and |f/>( 2 )| 2 , hereafter. Similar discussions are also possible for the negative incidence direction. The ^(s,qi,q 2 ) in 
the junction part is derived as 

*(A q u q 2 ) = 1 ]T(3/4r 1 (x 2 A + X 2 2 + 2^ cos ( 6 + ^( Ql - q 2 )) X hl X h2 ) , (53) 
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and 



1 X — > . , ,_• i / .) o £ / 27T 

\ A\^~ 1 I I "V^2 



i=l 



g x , ft) = - 2J3/4) 4 - 1 I + X/ >2 - 2- cos \6 - —(qi - q 2 ) 1 X M X, 2 1 , (54) 



where Xjj = Xi^(kRt, kr) is defined by eq.(|50|), with the circumference i?t of the tube into which the electron is 
transmitted. That is to say, Rt = R7 and kR t < kr < kR$ < 2tt for the positive incidence direction, while Rt = R5 
and fci?7 < kr < kR t < 2ir for the negative incidence direction. Note that the squared wave function \& in eqs.([5^) 
and (HJ) does depend on the sign of energy unlike the conductance in eq. (HJ) . Note also that the sign of E is inverse 
of that of e, E/\E\ = —e/\e\ = ±1, because the transfer integral 7 is negative. In the right hand sides of eqs.([5^) 
and (|54|), only the third term depends on the angle 9, and it becomes zero as E approaches zero, because X^i — ► 0. 
Therefore the spatial oscillation of 'J along the 9 direction decreases with decreasing \E\. 

There are several possible choice of the basic translation vectors, e\ and e 2 . The wave function represented by 
eqs.(|5^) and (|54|) should be invariant under the change of the choice. When e'i = e 2 and e'2 = e*2 — e 1 are used 
instead of (e'i, e'2), the corresponding x', y' axes are rotated ones by 7r/3 with respect to original x,y axes as shown 
in Fig.[l| Correspondingly, the polar coordinates and the labels are transformed as (r, 9) — > (r, 9 — 7r/3), (A, qi, q 2 ) — > 
(B, qi + q 2 , —qi — 1) and (B, qi,q 2 ) — > (A, q± + q 2 + 1, —qi — 1). For example, the site indicated by the square symbol 
near the origin O in Fig. [j] has the two different labels, which are (A, 0,0) defined by (ei,e 2 ) and (5,0, — 1) defined 



by (e'i, e'2). Under this transformation of (r, 9) and (s, q%, q 2 ), the values of eqs.(p3[) and (54) are invariant, that is to 
say, the results are independent of the way of the labeling, or choice of the vectors, e\ and e 2 . 

In order to relate the wave function to the transmission rate, the integral of the squared wave function defined 
below should be introduced: 

9 r p ,22 

$(kR t ,kr) = — d9(^(A- 1 r,9) + ^(B- 1 r,9)) = -J2J2^/^ 1X U kR ^ kr ) ■ ( 55 ) 
n ■'Q i=i j=i 

The third terms in r. h. s of eqs.(|53|) and (|54|) do not contribute to the integral because q\ — q 2 varies rapidly as a 
function of 9. We call 4> defined by eq.(|55|) a radial norm, hereafter. The radial norm is normalized so that it equals 
unity at r = Rt, i.e., at the exit of the transmitted wave from the junction part. It corresponds to the unit flow of 
the transmitted wave. Thus at the entrance of the incident wave into the junction part, the radial norm $ becomes 
the sum of the amplitude of incident wave 1/T and that of the reflected wave R/T = (1 — T)/T. Accordingly, the 
transmission rate T can be given by T = 2/($ + 1), where <!> = $(fci?7, k,R§) for the positive incidence direction 
and $ = &(kRs, kR-j) for the negative incidence direction. Figure ^| is the schematic development maps representing 
this relation between the radial norm and the transmission rate. It shows two junctions with a common thicker 
tube where the incidence direction is negative. Though they have different transmission rates, T\ and T 2 , their wave 
functions in the region CC'D'D and the transmitted wave are common. From now on, we call 2/($ + 1) an 'inverse' 
of the radial norm for simplicity. The inverse of the radial norm, 2/($(fei?7, kRs) + 1) ,2/(^(kR^,kRj) + 1) and 
the transmission rate T of eq.(pl]) are shown in Fig.|| as a function of E/E c by the dashed lines, the dotted lines 
and the solid lines, respectively. The curves are shown for the cases R7/R5 = 0.1,0.3,0.5,0.7 and 0.9. For each 
value of R7/R5, the inverse of the radial norm is close to the corresponding transmission rate. Therefore we can 
discuss the dependence of the transmission rate on R7/R5 and \E\/E C by the radial norm. For this purpose, we 
show in FigJ^ the radial norm for the negative incidence directions, ^(kR^, kr), as a function of r/R§ at six energies 
E/E c — kR 5 /(2n) — 0.1,0.3,0.5,0.65,0.8,1.0. There are four regions of the parameter space which have different 
characters; they are region I (1 > r/R 5 > 0.9; 1 > \E\/E C > 0) , II (0.9 > r/R 5 > 0.1; 0.7 > \E\/E C > 0), III 
(0.9 > r/R 5 > 0.1; 1 > \E\/E C > 0.7) and IV (0.1 > r/R 5 > 0; 1 > \E\/E C > 0). In the region I, all the curves in 
Fig.^ are very close to unity. It indicates that the corresponding transmission rate is near unity and independent of 
\E\/E C as is shown by the lines of R7/R5 = 0.9 in Fig|| In the region II, the radial norm shows monotonic decrease 
as either r/R$ or \E\/E C increases. The decrease with respect to r/R^ becomes steep as \E\/E C decreases and are 
almost proportional to (R$/r) 3 when \E\/E C < 0.1. In the region III, however, the radial norm oscillates as r/i?5 
increases. When \E\/E C = 1 the radial norm has a minimum value around r/R§ = 0.5. Correspondingly, the curves 
of i?7/i?5 = 0.5 are larger than those of R7/R5 = 0.3 and 0.7 at \E\/E C = 1 in Fig.^. It oscillates also with increasing 
\E\/E C when 0.4 > R§ > 0.2. The lines of R5/R7 = 0.3 in Fig.|^ shows the corresponding oscillations. Around the 
boundary between the region II and III, the decrease of the radial norm with increasing r/R§ becomes critical as is 
indicated by the plateau of the curve corresponding to \E\/E C = 0.65 in Fig.^. Correspondingly, all the curves for 
R 7 /R 5 = 0.7, 0.5 and 0.3 in Fig.| cross each other around the boundary energy E/E c = 0.7. In the region IV, the 
radial norm is nearly proportional to (i?5/r) 3 so that the transmission rate is almost proportional to (R7/R5) 3 at 
each value of \E\/E C . 







The oscillating behavior in the region III can be interpreted as a resonant effect in the following discussion. The 
wave function propagates both in the radial and the angular direction. Its wave length along the radial direction is 
assumed to be close to that of the plane wave obtained from the dispersion relation (|25|); the wave length is given 
by 27r/fc = R^\E C /E\. When the length of the junction part along the radial direction, i.e., R5 — Rj, coincides with 
the half of the wave length, the resonance occurs. From this discussion, the condition of the resonance can be easily 
obtained as 

R 7 /R 5 = 1-\E C /2E\ . (56) 

The minimum points of the curves in FigJ^ in the region III correspond to this resonance following approximately 
condition (|56|); it comes around r/R$ = 0.5 when \E\/E C = 1 and moves towards the left as \E\/E C decreases. Around 
the resonance point, the transmission is almost perfect. Such a resonance, however, does not occur in the region II 
and IV, since the radial norm becomes divergent when kr approaches zero, i.e., when either \E\/E C = kR^/(2n) or 
r/i?5 approaches zero. This divergence comes from the terms of krN%(kr) in eq.([55|) which is almost proportional to 
l/(fcr) 3 for small kr values. Accordingly, the transmission rate decays with the power law in proportion to (R7/R5) 3 
in the region II and IV. 



IV. COMPARISON BETWEEN THE ANALYTICAL RESULTS AND THE NUMERICAL RESULTS 



Since the effective mass theory is an approximation of the tight binding model as is explained in the section III, its 
results have to be confirmed by comparing them to the corresponding numerical results from the tight binding model. 
In this section, we shall compare the wave function calculated numerically by the tight binding model, Cj, with those 



determined analytically by eqs.(|53[) and (54), "J. For the numerical calculation with the tight binding model, the 
conditioned transfer matrix method explained in section II is used. Figure || shows the ratio of the corresponding 
quantities, {\cf}\ 2 + \cf} | 2 )/(T*) as a function of radial distance r, for the energy E = — 0.05|7|. The values are 
plotted for the (3,3)-(13,13) junction in Fig^-(a) and for the (10,10)-(13,13) junction in Fig^-(b). Here the incidence 
direction is positive; c^J is caused by the incident waves For the opposite incidence direction, similar results are 
also obtained. Here the factor 1/T = \h + \ 2 defined by eq.([49|) is necessary because the transmitted waves have unit 
flow in eqs.(p3|) and (Ejj), while the incident waves have unit flow in eq.(p^). One can see that the ratio is close to 
unit, i.e., the coincidence is fairly good. Correspondingly, the transmission rate by the conditioned transfer matrix 
and that by eq.([5l|) coincide well with each other when \E\ < 0.2|7|, as shown in Fig.||. The deviation from unity 
becomes smaller as \E\ approaches zero, because all the assumptions postulated in deriving eqs.(|5^) and ([54j) arc 
more valid for smaller \E\. The deviations, however, become larger when r ~ R5 in Fig.^-(a), while such a feature 
is not found in Fig.||-(b). The most natural explanation for this difference is that the deviations are caused mainly 
by the evanescent waves in the tube parts, and their effects become more significant as their decay lengths along 
the tube axis become larger, i.e., as they approach the extended states. The decay length is determined by ^7 in 



eq.(|27]) as aj ^ {2irl^ja/ R5.7) 2 — e 2 . Thus the discussion below is concentrated to the most extended evanescent waves 
corresponding to l 5 j = ±1. When i? 5 ^> R 7 as the case of FigJ|-(a), the evanescent waves in the thicker tube side 
are much more extended and have larger effects than those in the thinner tube side. As a result, the deviation near 
the thicker tube side is enhanced compared to that in the other place. On the other hand, when R§ ~ R 7 as the case 
of Fig.^-(b), the evanescent waves have similar decay lengths both in the thicker and in the thinner tube side, so that 
the difference of the deviation does not depend on the sides. 

To see the spatial variation of the wave function, we show \cj\ 2 of the (3,3)-(13,13) junction calculated by the 
conditioned transfer matrix method and the corresponding quantities determined analytically by eqs.(|5^) and ([54]). 
The former is shown by the closed symbols joined by the dotted lines and the latter is shown by the open symbols 
joined by the solid lines for E = d=0.2 |-y | in Fig.^-(a) and for E — ±0.23|'y] in Fig.||-(b). The horizontal axis is the 
number of the bar j, which are almost proportional to the radial distance r. The diamonds and the squares correspond 
to the positive energy and the negative energy, respectively. FigJ^-(a) shows that the results by the two methods agree 
quite well. Enhancement of the deviation at j = 16 ~ 20 is caused by the evanescent waves. In Fig.^-(a), the wave 
function does not show monotonic decrease along the incidence direction but has a peak structure. It gives rise to 
the resonant peak structures near E = =b0.2 j'y ] in Fig. |] which is already discussed in the previous section. In contrast 
to the good agreement in Fig.^|-(a), the deviation becomes quite large in Figj^-(b) especially near r = R§- 7 the wave 
function by the conditioned transfer matrix grows more rapidly than that by the analytical results with approaching 
the thicker tube side. When \E\ ~ E c , as in the case of Figjp-(b), the decay lengths of the evanescent waves are very 
large, and cause the large discrepancies. Corresponding to it, large discrepancies occur also in the transmission rate; 
only the data by the tight binding model show the sharp dips near E = ±0.237 in Fig.0. The sharp dip appears 
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approximately when 0.9 < \E/E C \ < 1, i.e., when the decay length of the evanescent waves in the thicker tube are 
larger than the diameter of the thicker tube. Evanescent waves with their decay length larger than the diameter of 
the corresponding tube are called 'quasi-extended' evanescent waves hereafter. 

In Fig|| data corresponding to sites belonging to the pentagonal defect or the heptagonal defect are shown by closed 
diamonds. They indicate that the numerical wave function calculated by the tight binding model is more localized 
at the defects than the analytical wave function. Hereafter a localization strength for the 8 direction at the defect 
is defined as the ratio between the numerical norm per site at the defect and that of the bars the defect belongs to. 
It is calculated by the tight binding model and represented by (J2i,keCj 

. \coA 2 + Kfc| 2 )(do + di)/(7|c | 2 + 7|c!| 2 ) at 
the heptagon and keCr \ci ly i\ 2 + \ci 1+ i y k\ 2 )(di 1 + di 1+ i)/(5\ci 1 \ 2 + 5|q 1+ i| 2 ) at the pentagon with the notation of 
section II. Here i £ C n means that site i belongs to the n membered ring defect. Since the radial distance r in the bar 
is almost constant, the analytical quantity corresponding to the localization strength has a maximum represented by 
ELiG^mI + l^,2|)7£i=i(l^,i| 2 + l^, 2 | 2 ) < 2 which is shown from eqs.(||) and (fj). Therefore we can say that 
the wave function is 'quasi-localized' at the defect when the localization strength is larger than two. The localization 
strength at the pentagon and that at the heptagon are shown as a function of the energy in Fig.|l0[-(a) and Fig.|lO|-(b), 
respectively. The solid lines and the dashed lines show the localization strength for the the (10, 10)-(13, 13) junction 
and that for the (3, 3)-(13, 13) junction, respectively. The solid lines shows that the the wave function is quasi-localized 
at both the defects when \E\ ~ E c and R7/R5 ~ 1, i.e., when there are quasi-evanescent waves in both the tubes. 
Furthermore they also show asymmetry with respect to the energy axis; the quasi-localization at the pentagon (at the 
heptagon) around E = —E c (E = +E C ) is larger than that around the energy with the opposite sign. These results 
indicate that the quasi-localization occurs due to the mixing between the quasi-extended evanescent waves and the 
'defect levels'. Here the 'defect levels' are defined as n discrete energy levels, — 2\<y\ cos(27r//n) (I = 1, 2, • • • , n), of an 
isolated n membered ring calculated by the tight binding model£3 Since the energy region around zero is considered 
now, the discrete level closest to E — 0, which is at E ~ — O.6I8I7I for the pentagon, and at E ~ +0.445|7| for the 
heptagon, is the most important. The closer distance in the energy between the two states enhances their mixing so 
that it causes the asymmetry with respect to the energy. Since the defect levels are caused by the discreteness of the 
lattice, these results can not be reproduced by the effective mass equation. 

The energy where the quasi-localized evanescent waves in the thinner tube side appear is around E = ±E c R^/Rt. 
When i?7/i?5 <C 1, it is distant from the considered energy region \E\ < E c so that there is no quasi-extended 
evanescent waves in the thinner tube side. Though there are the quasi-extended evanescent waves in the thicker 
tube side, the defect level of the heptagon cannot be mixed with them due to the large spatial distance between the 
heptagon and the thicker tube. This is the reason why the dashed line in Fig.|lCi|-(b) do not show the quasi-localization. 
In summary, we can say that the defect level is mixed with the quasi-extended evanescent waves only when they are 
close to each other both in the energy and in the space. 

In this way, the analytical results by the effective mass approximation are not appropriate when the evanescent 
waves become quasi-extended. The width of this energy region, however, is only about 0.1E C . Except these narrow 
energy region, the analytical results are appropriate enough. 

V. SUMMARY AND CONCLUSION 

In this paper, the junctions connecting the two metallic nanotubes with different circumferences by a pair of 
the pentagonal and the heptagonal defect are investigated. Both the wave function and the transmission rate are 
analytically obtained by the effective mass equation, and a close relation between them is found. To discuss the close 
relation, the junction part is divided into pieces according to the radial distance r from the thinner tube side and the 
radial norm of the wave function is evaluated in each piece. The radial norm is determined by the two parameters, 
\E\/E C and r/Rs where E c and R$ are the threshold energy and the circumference of the thicker tube, respectively. 
The transmission rate approximates to the inverse of the radial norm where r is substituted with the circumference 
of the thinner tube, R7. From the dependence of radial norm on the two parameters, the parameter space is roughly 
classified into the four regions; region I (1 > r/R 5 > 0.9; 1 > \E\/E C > 0) , II (0.9 > r/R 5 > 0.1; 0.7 > \E\/E C > 0), 
III (0.9 > r/R 5 > 0.1; 1 > \E\/E C > 0.7) and IV (0.1 > r/R 5 > 0; 1 > \E\/E C > 0). In the region I, the radial norm 
is close to unity so that the corresponding transmission is almost perfect, independent of \E\/E C . In the region III, 
the radial norm oscillates near unity as the function of either of the two parameters. The period of the oscillation 
with respect to r is approximately the same as half wave length. Here the wave length is that of the plane wave, i.e., 
it is obtained from the liner dispersion relation of the monolayer graphite, eq. (p5|) . Almost perfect transmission due 
to the resonance occurs in region III when the radial length of the junction part, R$ — R*?, coincides with this period, 
because the corresponding radial norm has a minimum value there. As either \E\/E C = kRc,/(2ir) or r/i?g decreases, 
however, the radial norm becomes divergent owing to the term krN 2 (kr) so that the corresponding transmission rate 
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approaches zero. It leads to the power Ja.w_decay of the transmission rate proportional to (R7/R5) 3 in the region II 
and IV instead of the resonant feature.uou 

The nanotube junction can be considered as a combined system, that is the quasi two-dimensional structure (= 
the junction part) connecting the two quasi one-dimensional structures (= the tube parts). The wave functions in 
the former part become the 2D waves, whose radial parts are given by the sum of Ji(fcr) and Ni(kr) (i = 1,2), while 
they are the ID plane waves in the tube parts. Nevertheless, the condition of the resonance in the region III, eq.([56l), 
is obtained from the wave length of the ID plane wave. It indicates that the four components of the wave function 
in the junction part, J^(fcr) and iVj(fcr) (i — 1,2), are combined appropriately to be fitted well with the ID plane 
waves in the tube parts. On the other hand, the wave function shows the power law decay along the radial direction 
in the region II and IV, because only the component N2(kr) becomes dominant. In other words, the wave function 
in the junction part has the two dimensional character in the region II and IV so that its matching with the plane 
wave in the tube parts becomes worse than in the region III. The wave function in the junction part has different 
dimensionalities in this way. 

The wave functions analytically obtained by the effective mass equation are compared with that numerically obtained 
by the tight binding model. When \E\ is very close to E c , the discrepancies between them becomes considerable. From 
the dependence of the discrepancies on the radial distance, we speculate that they are mainly due to the evanescent 
waves with their decay length larger than the diameter of the corresponding tube and they are called quasi-extended 
evanescent waves. The quasi-extended evanescent waves cause the sharp dip of the transmission rate, which is absent 
in the analytical result. When the quasi-extended evanescent waves and the energy levels of the n membered ring 
defect (n — 5, 7) are close both in the energy and in the space, mixing between them appears as the enhancement of 
the wave function at the corresponding defect. Nevertheless coincidence between the analytical wave functions and 
the numerical ones are fairly good except the narrow energy region 0.9 < \E/E C \ < 1. 

The nature of the wave function discussed in this paper can be observed by Scanning Tunneling Microscopy (STM) , 
because STM images reflect the Local Density of States (LDOS), which are proportional to the sum of the squared 
wave functions, eqs. ( |53| ) and (|54|), in the nanotube junction part. The third terms in eqs. ( |53| ) and (|54|) give rise to the 
y/3x -s/3 pattern or the oscillation along the 9 direction in the STM images. They fade away as \E\ approaches zero. It 
is also expected that the images depend on the direction of the electronic current 'along the tube axes', which should 
not be confused with that of the tunneling current 'from the STM tip to the sample'. When Rs/R? ^S> 1 and E ~ 0, 
for example, the squared absolute value of the wave function in the junction part is almost proportional to r 2 when 
the current flows from the thicker tube to the thinner tube and to 1/r 4 when it flows along the inverse direction. We 
expect that these results will promote the investigation of the nanostructures composed by the nanotubes including 
the nanotube junctions. 
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APPENDIX 

In the tight binding model, the probability flow from site i to site j is represented by Im(a*Hijaj), where is 
the amplitude of the wave function at site i and Hij is the hopping integral connecting site i and j. Conservation 
of the flow is guaranteed by the tight binding equation, Eai = JV Hijaj. In this paper, ifjj is nonzero only 
when i and j are nearest neighbors. Consider bars for j < forming to the (to, n) nanotube part in Fig.0. The 
amplitudes of wave functions are represented by r^-.i and in the same way as Fig.^[ Then the flow between 
j — l'th bar and j'th bar is represented by 7 Eti" I m ( r j-i Hereafter the common hopping integral between 

the nearest neighbors 7 and the lattice constant a is chosen to be units. The flow J corresponding to ipn+ m ( |26| ) is 
J = TOsin(?7 — &2 + 27r/3) + nsin(— 77 + k± + 27r/3), where the first term comes from the m bonds, 1 < i < m, and 
the second term comes from the other n bonds. Here ki — k ■ e*j and 77 is the angle of Rj with respect to x axis, i.e., 
(cos(?7), sin(?7)) = (\/3(to + n)/(2Rj), (n — m)/(2Rr)) and R7 = \J m 2 + n 2 + ran. By developing about ki to the first 
order and using the boundary condition mk\ + nk2 = 0, the flow J is represented as 

J = Rr - ^-rn(n + m) + 0(k 2 ) + 0{k 2 2 ) . (57) 
2R7 
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The absolute value of the second term in ( j57| ) is less than \/3ttRt/ R$, since the energy region where the channel 
number is kept to two, i.e., kR$ < 2?t, is considered now. It can be deduced from it that the flow is the almost 
constant value, R-?, when R 5 ^> y/Sir. On the other hands, the value of |a7+| 2 obtained by substituting F's in eq.Q) 
with those of £q.(26) for ip4£+ is also Rr. Therefore the normalization factor in eq.(41) is correct. The normalization 
factor of Rcf.E3 and Ref.Ea is different from it by a factor \/2. This difference does not change the transmission 
rate, but when ratio between the wave function in the tube parts and that in the junction part is considered, the 
normalization factor in eq.([i"l"|) have to be used. 
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FIG. 1. The development map showing the bond network of the (2,2)-(2,5) junction. The filled circle at the bottom and 
the open circle at the top in each bar indicate an identical atom. To form the junction, the development map is rolled up so 
that the filled circle and the open circle coincide with each other in each bar. The solid lines and the dashed lines represent the 
bonds within the bar and those connecting the neighboring bars, respectively. The points Ps(= Qs), Pt{= Qi) and O are the 
centers of the pentagonal defect, that of the heptagonal defects and the origin of the coordinate, respectively, which are shown 
also in Fig. ^| The two sets of the translation vectors, {ei,e*2} and {e'i,e'2}, and corresponding sets of orthogonal vectors 
{e x , e y } and {e x i , e y /} are shown. Note that definition of sublattices A and B are changed between the two sets. 

FIG. 2. Notation to represent the amplitudes of the wave function. It is illustrated with the fist bar and the second bar in 
Fig.jl]. The site in the j'th bar can be classified to the two groups according to whether it connects with the right (j + l)'th 
bar or with the left (j — l)'th bar. The former and the latter are denoted by rj,i (closed triangles) and lj t i (closed squares), 
respectively, when the corresponding site has the bond i, which connects the neighboring bars and is numbered from the bottom 
to the top. 

FIG. 3. Development map of the nanotube junctions. The lines 'E7P7P5E5 ' are connected and become identical with the lines 
'F7Q7Q5F5', respectively. The rectangles 'E7P7Q7F7' and 'P5E5F5Q5' form the thinner tube and the thicker tube, respectively. 
The triangle 'OP7P5' is the same as the 'OQ7Q5' rotated by 60 degrees. The quadrilateral 'P7P5Q5Q7' forms a junction part 
with a shape of a part of a cone. A heptagonal defect and a pentagonal defect are introduced at Pt(= Q7) and Ps(= Qs), 
respectively. The direction of the circumferences of the tubes in the development map is represented by their angles 775 and r\i 
measured anti-clockwise with respect to the x axis which is defined in Fig. [j]. The angle between the axes of the two tubes, (f>, 
is defined as <f> = 777 — rj 5 . 

FIG. 4. Transmission rates per channel, T, as a function of the energy. The transmission rate calculated by eq.(|5l|) are 
shown by the solid lines, and those calculated by the tight binding model are shown by the dotted line for the (10,10)-(13,13) 
junction and by the dashed line for the (3,3)-(13,13) junction. The considered energy region is \E\ < E c where the channel 
number is kept to two in both the tube parts. The threshold energy E c for the effective mass theory is about 0.242|'y| and that 
for the tight binding model is about 0.239|7|, where 7 < is the hopping integral. One can easily see the agreement between 
the two methods is fairly good. 

FIG. 5. Schematic development maps representing the relation between the radial norm "I> and the transmission rate T. 
The thicker tube and the region CC'D'D are common in both the development maps 

FIG. 6. Inverse of the radial norm for the positive incidence direction defined by 2/($(fei?7, kR$) + 1) , and that for the 
negative incidence direction 2/(&(hRs, kRr) + 1) are shown by the dashed lines and the dotted lines, respectively. Here the 
radial norm $ is defined in eq.(p5|). Values of R7/R5 are 0.1, 0.3, 0.5, 0.7, and 0.9, which are attached to the corresponding 
lines. The transmission rate per channel calculated by eq.(|5l]) is shown by the solid lines. It can be seen that the radial norm 
gives a good estimated value of the transmission rate. The horizontal axis is the energy normalized by E c . Here \E\ < E c 
represents the energy region where the channel number is kept to two in both the tube parts. 

FIG. 7. Radial norm for the negative incidence direction, i.e., <&(kR$,kr) is shown as a function of r/R$ for five energies, 
kR 5 /(2-K) = E/E c = 0.1,0.3,0.5,0.65,0.8 and 1.0. Here $ is defined by eq.(||). The values of E/E c are attached to the 
corresponding lines. Both the axes are represented with log scale. Inset shows the regions of the parameter space, I, II, III and 
IV, which are used in the discussion in the text. 

FIG. 8. Comparison between the wave functions \c^\ 2 calculated by the tight binding model, , and ^ defined by eqs.(^) 
and (pi]); It shows flc^ | 2 + \c^ 2 J | 2 )/(T*I/) as a function of radial distance, r, (a) for the (3,3)-(13,13) junction and (b) for the 
(10,10)-(13,13) junction. They are shown for j = 1, 2, • • • , h. The energy is E — — 0.05|'y| . Here \c^i\ 2 is caused by the incident 
waves with positive direction, i/>jf + . Here 1/T = \h+\ 2 defined by eq.(^) is used for the normalization. The vertical axis is 
represented with log scale. Data at the pentagon and at the heptagon are shown by the closed diamonds. 
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FIG. 9. Norm of the wave function of each bar, \cj\ 2 , of the (3,3)-(13,13) junction (a) for E = ±0.2|7| and (b) for 
E = ±0.23|7|. The horizontal axis represents the position of the bar j shown in Fig.l. The incidence direction is positive, i.e., 
from the right to the left in this figure. The dotted lines with closed symbols and the solid lines with open symbols correspond 
to those calculated by eq.(|l3|) and those calculated by eqs.(^) and (p^), respectively. The diamonds and the squares correspond 
to the positive energy and the negative energy, respectively. 

FIG. 10. The localization strength for the 9 direction (a) at the pentagon and (b) at the heptagon as a function of the 
energy. The incidence direction is positive. The solid lines and the dashed lines correspond to the (10, 10)-(13, 13) junction and 
the (3, 3)-(13, 13) junction, respectively. 
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